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1. INTRODUCTION 

Human demands for energy will not decrease so we are going to run out of fossil fuels for energy. 
Nuclear energy is not likely to be a major source of world energy consumption because of the relative risks 
associated with unleashing the power of the atom and we have no choice but to invest heavily in renewable 
energy production. Solar energy is one of the most promising renewable sources of electricity power in 
the world. Moreover, the coordinated inter-connection of solar energy sources between global electricity 
markets can supply more flexibility and balancing to the grid. This is why the grid-connected photovoltaic 
systems have a great bright future. 

Achieving the reduction of PV system manufacturing costs remains a major task. Among the trends 
to achieve such a reduction is the elimination of the DC/DC converter stage. Conventionally, the first 
converter stage, which is usually placed between the PV arrays and the inverter [1]-[3], achieves the MPPT 
whereas the inverter stage delivers and controls the energy injected into the grid. Therefore, to achieve this 
cost reduction, the three-phase inverter must also take care of the maximum power point tracking. 

Maximum power point tracking is mandatory to maximize photovoltaic systems efficiency. To this 
end, several MPPT control strategies have been largely published in the last few years [1]-[8]. Perturb and 
observe (P&O) and incremental conductance (IC) algorithms are the most widely proposed in literature. 
They are the simplest algorithms to implement [3]-[5], but the dilemma -the choice between convergence 
speed and output fluctuations [5], [6], has led in recent years to several research aimed at improving these 
two techniques [9], [10]. The nonlinearity of photovoltaic characteristics makes the control of PV systems by 
conventional control strategies a complex task. However, the recent involvement of nonlinear controls has 
enriched the field of research and has proved most appropriate for the control of nonlinear systems [7], [8]. 
In particular, [6] and [11] are comparative studies between the backstepping control and other conventional 
controls which clearly showed the performance of the backstepping control. 
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In this paper, two controllers for the three-phase inverter, based on backstepping approach, are 
proposed and discussed. The first control objective is to track the maximum power point (CPp/OVP=0). 
The second objective is to ensure a grid connection with unity power factor (output current must be in phase 
with the grid voltage). The proposed controllers are developed in the synchronous orthogonal frame and 
the voltage phase angle of the grid utility is detected using a three-phase PLL (Phase-Locked Loop). 

This paper is organized as follows. In section two, the system description and the mathematical 
model are presented. Section three develops the control strategies of proposed system. Simulation results are 
presented and compared in section four. In the last section, a conclusion followed by the reference list. 


2. DESCRIPTIONS AND MODELING OF PV SYSTEM 

It is essential to find a simple cell model that would be useful for PV application. The PV cell is 
basically a p-n semiconductor junction that converts light energy into electrical energy. The PV cell can be 
represented by a current source controlled by voltage, sensitive to temperature and solar radiation, in parallel 
with one diode and a shunt resistor Rs, and the whole is in series with a resistor R, (Figure 1). 
The mathematical model of solar cells is detailed in [3]. Although it is a simplified model, this equivalent 
circuit is sufficiently accurate to represent the different types of photovoltaic cells. According to [12], both 
resistors (Rs and Rsh) may be neglected. However, it is demonstrated that the series resistor (Rs) has 
a significant impact on the accuracy, between the maximum power operating point and the open circuit 
voltage, of the current-voltage characteristic curve inclination. 

The photovoltaic array is a multiple associations, in series and parallel, of PV cells. The PV 
generator considered in this paper is composed by thirty-three SM55 Siemens panels connected in series. 
The electrical specifications for one panel are enlisted in Table 1. The modeling of the SMS55 panel, using 
Matlab/Simulink, allowed tracing its characteristics for different values of irradiance and temperature which 
are shown in Figure 2 and Figure 3. Figure 4 shows the power-voltage characteristics of the PV generator 
considered in this paper under changing climatic conditions. The coordinates of the Maximum Power Points 
(MPPi) are summarized in Table 2, and will be used for the verification of the simulation results. 


Table 1. Electrical Specifications for SM55 Solar Panel 


Parameter Value 

Maximum power 55W 

Current at the maximum power point 3.15A 

Voltage at the maximum power point 17.4V 

Maximum current (short circuit output) 3.45A 

Maximum voltage (open circuit) 21.7V 
Current temperature coefficient 1.2mA/°C 

Number of series cells Ns 36 
Number of parallel modules Np 1 
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Figure 1. Equivalent circuit of PV cell Figure 2. Power-voltage characteristics for the SM55 
panel at 1000W/m2 with varying temperature levels. 
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Figure 3. Power-voltage characteristics for the 


SM55 panel at 25°C with varying radiation levels. Generator 


Table 2. Maximum Power Points (MPP) In Figure 4 


Maximum power point Voltage[ V] | Power [ W | 
MPP1 584 1847 
MPP2 531.6 1669 
MPP3 573.2 1447.8 


The basic power circuit of proposed single stage PV-system (Figure 5) consists of a DC-link 
capacitor which is directly connected to the three-phase IGBT inverter. Grid connection was performed 
through a low-pass filter used to reduce the ripple components due to the switching actions in PWM inverter. 
It is assumed that all three phase inverter switches are ideal, the low-pass filter phases are identical 
(Li=L2=L3=L & Ri=R2=R3=R) and the grid voltage is symmetric. The three-phase model of the inverter with 
its low-pass filter is detailed in [7]. A simplified model can be obtained in the synchronous orthogonal frame 
rotating at the angular frequency of the grid voltage. For this purpose, the power-invariant qd-transformation, 
from balanced three phase electrical quantities to balanced two phase quadrature quantities, has been used. 
The abc-qd transformation matrix is given below (1): 


cos(@) cos(@— =) cos(@ + 5 


T@) = |= 1 
K Ji as —sin(é — =) —sin(@ + =) S 


where @ is the angular position of the dq-frame. The state-space model can be re-written in the new 
reference frame, as (2) [7]: 


Ls = Ig + ol- + Ky 


é —R E V 
le =at 1g + |, (2) 


Vp = (Ip — laKa — [qKq) 


where:( I»; V, ): PV array current and voltage;(Ka; K,): Direct and quadrature inverter control inputs; 
(Ia; I,): Direct and quadrature injected current components; (Ea; Eg): Grid direct and quadrature voltages; 


[Ka Ko] =T@)L kt ke ks f; [la K| =TOLi iz is I": [Ea Ey] =T@Ler e, ez J") 


3. CONTROLLER DESIGN 

In order to match with the proposed modeling presented earlier (2), the three-phase grid currents and 
voltages are transformed into direct and quadrature axis components. The controller outputs (Kg & K, ) are 
then transformed into three-phase components using the inverse dq-abc transformation, and then a PWM 
control is used to make them suitable for switching the inverter switches (Figure 6). 

The injected currents have to be synchronized with the grid voltages. To this end, the grid voltage 
phase angle is detected using a 2™ order phase locked loop (PLL). The structure of the PLL implemented in 
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this work (Figure 7) uses the grid voltage abc-dq transformation to track the grid voltages phase angle. A PI 
regulator was used to generate a corrective phase angle (6st), that is fed back into the grid voltage abc-dq 
transformation module, from the quadrature voltage error [13]. The block diagram representation of 
the backstepping control approach is illustrated in the following Figure. 
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Figure 5. PV inverter with the low-pass filter 
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Figure 7. Basic principle of PLL 


In the synchronous d-q reference frame (£,=0), the injected powers are decoupled and simplified (3): 


| P = Eala (3) 


P = Eala + Eq! 
f Q = —Ealy 


Q = Egla — Eql 


According to (3), active and reactive powers can be controlled by direct and quadrature current 
components (I4 & I4) respectively. In order to achieve a minimum injection of reactive power the quadrature 
current reference Igref must be set to zero. If we neglect ohmic loss, the power conservation principle 
gives (4): 
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N L AR 

laret = a Ea (4) 
Finally, when the PV-generator power is at its maximum state, (Figure 3), its derivative with respect 

to PV-voltage is zero. Table 3 summarizes the selected dynamic outputs and their references for each 

controller. 


Table 3. Dynamic Outputs and their References 














The dynamic output considered Output Reference 
First Controller OP, e 
yı Ty ay ref 0 
y=Bil=| a m| IY JLo 
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Second Controller 7 OP, Py 
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where: eq = Iq — lareg 18 the direct current error, eg = Iq — Igrer 18 the quadrature current error. 

The output chosen for the second controller is simpler than the first, but it is a challenging and 
time-consuming task to tune the appropriate backstepping parameters in this controller. This is largely due to 
the fact that establishing stability for switched systems 1s difficult. Therefore, in the first controller, we have 
chosen to also involve the direct current component 74 in the control of the output power. 


3.1. Formulation of MPPT control law (controller 1&2) 
Let us define the first tracking error €, and its LFC (Lyapunov Function Candidate) V}, using (2), it 
is possible to deduce (5-6): 








al ð?I ôl. f (5) 
A = ae a p p = 
E€ = V1 Yıref = lp + bay => a = (Mara tge) = E -laka — hka) 
1 . . Ef (6) 
Vie e = Vy = €,€, = =U — laKa —1,Kq) 
2 
where: f = V, r F2 ze . The stabilizing controls Kg and Kg is chosen as follows (7): 
p p 
Bı + AKa + A2K, = 0 (7) 
_ lof À — laf, _ qf ; es 
where: B, = = + ¢C,€,;A,= Pa A, = = and c4 is a positive design parameter. 
With the above choice V, = —c,e? becomes defined negative, €, is proved to be asymptotically 
stable and converge to zero by the Lyapunov design. This means that y4 —> Yıref , so P, — (P,) upp- 
3.2. Formulation of output power control law 
3.2.1. First controller 
According to Table 3 we define the error €, and its derivative are deduce from (2) will be (8): 
E2 = Y2 — Yzref 
: -R Eq V . -R V 
la = 2e, (Zhu + wl, -+2 Ka — larer ) teolt = hue KR) (8) 


: È . 
Where: lares = 5° = = ý= = (lp — laKa — GK.) 
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Let us consider the following new LFC and its derivative with respect to time deduce frome (8) 
will be (9): 


V, = Vi + e? (9) 
Vp ; -R Vp 
Vz = —cye? + 2ez lea i Pole = Ld 4 7 Ka ` laref) + eg} ola + lq + = Ka )] 


Then the stabilizing controls for the second step (Ką and Kg ) is chosen as follows (10): 


B> + A3Kq + A,K, = 0 (10) 
where: 

Vp E€ Vp E R 
A3 = ea 7 te la» Ayse ta atb t o€ + ea (— la + ol u. p) + eg (—wlag —=1q) €z 
is a positive parameter. With the above choice V, = —c,e% — c,eŻ < 0, eand €, are proved to be stable and 


converge to zero by the Lyapunov design. This means that y converge to y,.~-. We are finally in a position to 
determine the stabilizing control signals, Kg & Kq, from (7) and (10) will be (11): 


A, | a Pa a 4 —_ A, ai E (11) 
A, A4] lKa B, => K| [4 4,] [Bo 


3.2.2. Second controller 
The tracking error €, 1s reduced to (12): 


From (2) ; —R E Vp 
E2 = lq- lge =q => é, = —Wlg +1, -1+ *K, (12) 


And the control low K, can be extracted directly from the second equation of mathematical 
model (2), by choosing Kg such that (13): 


Ky = (Cal Folk = 14) © €&=-M6 DB V, = V, +66, = —c 6? — ce? (13) 


Using (7) and (13) will be (14): 


ee-ie S] [aout 1 
EN =I, 


4. SIMULATION RESULTS 

The theoretical performances of the proposed controllers discussed in section three will be 
illustrated by simulation in this section. The PV system (Figure 6) is simulated jointly with each controller, 
using the instantaneous three phase model, in Matlab/Simulink environment (Figure 8). The model in d-q 
axis (2) is only used in the controllers design. Important simulation parameters are given in the Table 4. 
Controllers’ parameters have been selected using a ‘trial-and-error’ search method. In order to prove 
the robustness of the control algorithms, the simulation is performed with the following scenario: 
-A temperature increase from T=25°C (298.15K) to T=45°C (318.15K) after Isec of start of simulation, then 
returns to T=25°C at 1.6sec, as shown in Figure 9. 
-A solar irradiation drop from 1000W/m? to 800W/m? after 0.4sec of start of simulation, then returns to 
1000W/m? at 0.3 s as shown in Figure 10. 
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Table 4.Characteristics Of Controlled PV System 


Parameter Symbol Value 
PV-array power Pp 1847W 
DC bus capacitor C 100u F 
Switching frequency Fs 10kHz 
AC source Vg 110V 
Line frequency f 50Hz 
Filter parameters L 12.5mH 
R 0.5655 Q 
controller parameters Cy 8x104 
C2 9,7x10* 
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Figure 8. Simulation schema for the proposed PV system 
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Figure 9. Temperature changes 
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Figure 10. Solar irradiation changes 


The first and second controllers are based respectively on (11) and (14), and use the same 
backstepping parameters (Table 4). Simulation results of the PV array power and voltage under transient 
condition are shown in Figure 11 and 12, which correspond very well to the MPP coordinates summarized in 
Table 2. 
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Figure 13 illustrates the change of direct and quadrature current components. It can be clearly seen 
that the second controller can track the reference values with a fast transient response, but the first controller 
is the most accurate and the least disruptive of the grid. 
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Figure 11. PV power and DC-bus voltage behaviors 
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Figure 12. Zoomed portion of Figure 11 
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Figure 13. Current components behavior 


Figure 14. shows the grid voltage and the injected current obtained with each one of the proposed 
controllers (the grid current scale was multiplied by 10). 


It can also be observed that the amplitude of the phase current changes when the solar radiation or 
the temperature are varied. It can easily be seen in the zoomed portions (Figure 15) that the grid current is 
sinusoidal and in phase with the grid voltage which proving that the power factor unit is well achieved. 
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Remark: It should be noted that it is very difficult to tune backstepping parameters of the second 
controller to establish the stability, unlike the first controller whose parameters influence precisely on 
the precision and settling time. Although the current of the first controller has fewer ripples compared to 
the second, the use of a second-order filter (LCL for example) will surely improve the quality of 
the injection. 
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Figure 15. Zoomed portions of Figure 14 
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5. CONCLUSION 

In this paper, two Backstepping controllers have been developed for a three phase single stage grid 
connected photovoltaic system, without using the conventional DC/DC converter for MPPT control, in order 
to ensure the operation of the PV system at maximum power point with unity power factor. Using both 
theoretical analysis and simulation, it has been proven that the proposed controllers, although they require 
more computation, have given accurate results with virtually zero errors, it has also been deduced that 
the correct selection of system outputs helps to facilitate control tasks. Future work will deal with 
the mitigation of switching noise using LCL filter 
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